Atomic number fluctuations in a mixture of two spinor condensates 
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We study particle number fluctuations in the quantum ground states of a mixture of two spin- 
1 atomic condensates when the interspecies spin-exchange couphng interaction C12/? is adjusted. 
The two spin-1 condensates forming the mixture are respectively ferromagnetic and polar in the 
absence of an external magnetic (B-) field. We categorize all possible ground states using the 
angular momentum algebra and compute their characteristic atom number fluctuations, focusing 
especially on the the AA phase (when ci2/3 > 0), where the ground state becomes fragmented and 
atomic number fluctuations exhibit drastically different features from a single stand alone spin-1 
polar condensate. Our results are further supported by numerical simulations of the full quantum 
many-body system. 
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I. INTRODUCTION 

Since the first production of an atomic ^^Na conden- 
sate in an optical trap [1| , spin degrees of freedom for con- 
densed atoms become accessible, which has since given 
rise to a rich variety of phenomena such as domain for- 
mations U , spin mixing dynamics Q , topological defects 
[J|, etc. The properties of a three-component {F — 1) 
spinor condensate are first studied by Ho i5|] and Ohmi 
Q . Many predictions are experimentally verified [31 , the 
most fundamental property concerns the existence of two 
different phases: the so-called polar and ferromagnetic 
states, respectively corresponding to the F = I state of 
^■^Na and ®^Rb atomic condensates. 

In contrast to a scalar condensate, both spatial and 
internal spin part of the wave functions are required to 
discuss a spinor condensate. For both F = 1 state of 
^■^Na and ^^Rb atoms, density-density interactions are 
significantly larger than the spin-exchange interactions. 
The single spatial-mode approximation (SMA) ^5-9] is 
often adopted, whereby one adopts mean field approxi- 
mation (MFT) to determine the condensate spatial wave 
function neglecting their spin dependence. The spin de- 
grees of freedom is then considered assuming the spa- 
tial wave function is identical. The spin related prop- 
erties can be investigated use either a "mean-field" the- 
ory or "semiclassical treatment" 0, Q , where the spinor 
is described by a vector formed by three c numbers 
Ca(0(" = 1,0,-1), or use many body theory 0, |1] 
treating the spinor from the bosonic mode as operators 

Qa, satisfying aa,al = 6aj3- The quantum ground 

states for a spin-1 condensate have been studied exten- 
sively. 

It was initially predicted that the ground state of 
^^Na BEC (c2 > 0) is either polar (no = N) or anti- 
ferromagnetic {ni = n_i = N/2) in the mean-field the- 
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ory [1, Q . A quantum treatment based on the SMA by 
Law, Pu and Bigelow [3,|M| revealed, however, the ground 
state for ^"^Na atoms is actually a spin singlet with prop- 
erties (ni = riQ = 7T._i = N/3), drastically different from 
polar state predicted within the mean field theory. Fur- 
ther studies pointed out that this spin singlet state is a 
fragmented condensate with anomalously larg e number 
fluctuations and thus has fragile stability [ifl, [ill • The 
remarkable nature of this fragmentation is characterized 
by three macroscopic eigenvalues (see above) for the sin- 
gle particle reduced density matrix, which is capable of 
exhibiting anomalously large atom number fluctuations 
Ani,o,-i-A^. 

The interests in spinor condensates extend to higher 
spins [12, [l3| and quantum ground states are already 
well known and categorized for spin- 2 condensates [14 |. 
Atomic Feshbach resonance was implemented in a double 
condensate, enabling tunable interactions, whose effects 
on superfluid dynamics and controlled phase separations 
are observed 15]. The spin exchange interaction between 
individual atoms can be precisely tuned through "optical 
Feshbach resonances" [lg| by adjusting the two s-wave 
scattering lengths ao and 02 . This inspired several recent 
theoretical studies on mixtures of spinor condensates |17l - 
|20 | and tunable or controlled dynamics |2l| . 

In this paper we report anomalous fluctuations for the 
numbers of atoms in a binary mixture of spin-1 conden- 
sates. We hope to stimulate experiments, using the most 
relevant experimental case, the mixture of ^"^Na (polar) 
and ^^Rb (ferromagnetic) condensates in their F — 1 
manifold, as an example. The quantum spin properties in 
the special ground state of the AA phase, where the inter- 
species anti-ferromagnetic spin-exchange is large enough 
to polarize both species but forming an maximally en- 
tangled state between two species, are studied and we 
give the exact number fluctuations distribution. Then we 
resort to numerical diagonalizations to show that parti- 
cle numbers and number fluctuations undergo dramatic 
changes as inter-species coupling ci2/3 varies. 




FIG. 1: (Color online) The dependence of ground-state phases 
on /3 at fixed values of ci/?i — —1, C2/32 — 2. All interaction 
parameters are in units of |ci/3i|. 



II. THE MODEL HAMILTONIAN FOR THE 
MIXTURE 

Intra-condcnsate atomic interaction takes the form 
V,(r) = (aj + (3jFj ■ Fj)S{r) with j = 1, 2 for the ferro- 
magnetic (/3i < 0) and polar (/32 > 0) atoms. The inter- 
species interaction between the ferromagnetic and polar 
atoms is described as Vi2(r) — ^{a + (3Fi •F2 -|-7Po)^(r) 
[17| , which is more complicated because collision can oc- 
cur in the total spin i^tot = 1 channel between different 
atoms, in contrast to intra-condensate interactions be- 
tween identical atoms [13,1241 ■ The parameters a, /3, and 
7 are related to the s-wave scattering lengths in the the 
total spin channels T^, analogous to spin-1 condensates 
[5|, |6| . Pq projects an inter-species pair into spin singlet 
state and ^ = AIiM^/iMi + M2) denotes the reduced 
mass for the pair of atoms, one each from the two differ- 
ent species with masses Mi and M2 respectively. 

Our model Hamiltonian is given by 

H = H1+H2+H12, (1) 



-yW-WjWjWi 



2 
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H12 = i /"dr{a*l$]$,*, 

+l3i>mFui ■ F2,fe$fc*, + |oto| . (2) 

H2 is identical to Hi except for the substitution of sub- 
script 1 by 2 and ^i by $i. The latter two are atomic 
field operators for the spin state |l,i). O = ^i$_i — 

We adopt the SMA [3-Q for each of the two spinor 
condensates with modes ^(r) and $(r), i.e., setting 

^i = Oi*, $, = 6,$, 

with Hi (bi) the annihilation operator for the ferromag- 



aj,a] 



netic (polar) atoms satisfying [ai,aj] — and 

6ij (and the same form of commutations for bi). The 
spin-dependent Hamiltonian for our mixture mode then 
reads 



H 



ClA-f,2 , C2^2-f,' 



H n — Fl ■ F2 



Cl2TAt 



6 



612612, 



(3) 



where Fi = aJFujaj (F2 = 6|F2ij6j) are defined in 
terms of the 3x3 spin-1 matrices Fiij (F2ij), and 

Q\, = albl - a\bli - alibi 

creates a singlet pair with one atom each from the two 
species, similar to the following two 



it = (4)^ 



2a\dli, 



?t _ fhU2 



tf,t 



B^ = ibir-2bibii 



for intra-species spin-singlet pairs [lO[. The interaction 
coefficients are ci = /dr|^(r)| , C2 = /dr |$(r)| and 
C12 = /dr|^(r)| |<I>(r)| , which can be tuned through 
the control of the trapping frequency. Here we focus 
on the spin-dependent part when we assume that the 
two species are sufficiently overlapped. The scattering 
properties between any pairs of atoms in specific Zee- 
man hyperfine component states are determined from 
their corresponding interaction potentials observing the 
symmetries of the two states. We notice that under the 
so-called degenerate internal-states approximation (DIA) 
the 7 term vanishes [23|. Within the DIA, the Zeeman 
hyperfine state of an alkali atom is expanded in terms 
of the electronic (valence electron) and its corresponding 
nuclear spin states. Between any two atoms, the cor- 
responding interaction potential is constructed from the 
appropriate weighted spin singlet /triplet potentials when 
expanded out correspondingly. When the spin singlet 
and triplet potentials are further approximated by low 
energy scattering pseudo-potentials, the total scattering 
properties are determined completely by two scattering 
lengths of the singlet and triplet potentials. The param- 
eter 7 is therefore no longer needed and the Hamiltonian 
([3]) finally reduces to 



H, 



£]Ap2 



£2^p2 



^Fi.F2. 



(4) 



III. THE GROUND STATES 

We rewrite the the Hamiltonian ^ as 
H = afj + bPl + cF'^ , 



(5) 



with a = ci/3i/2 - ci2/3/4, b = C2;92/2 - ci2/3/4, and 
c — C12/3/4, F = Fl -I- F2 is the total spin operator. 



The eigenstates of ([5]) are the common eigenstates for 
the commuting operators F^ , F^ , F^ , and Fz , given by 

\F,,F2,F,m}= J2 C^Cmv,F.^rnjPurn,}\F2,m2}, 

7ni7n2 

with the uncoupled basis states 

\Fi,m,)^Z;^F,-f^-"'^{a\f^iA^)^^^-^^y'\0), 

and anafogously for \F2,m2) whic h sp an a Hilbert space 
of dimension {Nj + l){Nj + 2)/2 ^. C is the Clebsch- 
Gordon coefficient, Zj is a normahzation constant, and 
Fj- is the fowering operator for rrij. The corresponding 
eigenenergy is 

E ^ aFi{Fi + 1) + bF2{F2 + 1) + cF{F + 1). (6) 



Given Nj , the aUowed values of F,- are Fj = 0,2,4, ■■■ Nj 



if Nj is even and Fj = 1, 3, 5, 
\Fi~F2\ ^F^Fi+F2. 



' Nj ii Nj is odd, satisfying 



We next consider the special case of TVi = A^2 = ^ and 
for N even. Minimizing the energy (|6]), we can get the 
ground state phases determined by different parameters 
a, b, and c, or equivalently Ci/3i, 02/^2 , and Ci2/3. Taking 
N = 100 for example. Figure [T] presents the results for 
the order parameters (Ff), (F^), and (Fi -^2) = ((F^) - 
(Fl) - (F^))/2 Irf] as C12/3 changes but fixed ci/3i = 
— 1, C2/32 = 2 (in units of |ci/?i|), which are found to 
agree very well with the mean field ones obtained from 
simulated annealing [17l |. 



Phases 


Parameter range 


Ground states jFi , F2 , F, F) 


FF 


-00 < C12/3 < Z1M^C2P2 


Cl''.^^^.^\N,N)\N,N) 


MM- 


'<T"c2/?2<Ci2/?<0 


C^^'-J:^^X'^N,N)\F,,F,) 


MM+ 


< C12/? < ^C2/32 


E C^,;.^^>^:™?|7V,mi)|F2,m2> 

mi ,7712 


AA 


^C2/32 < C12/3 < +00 


J:ct}%,;N,-m,\N,rm)\N,-rm} 

mi 



TABLE I: Quantum phases for the ground state. 



r 



Overall, there are four different phases in the 
general case, separated by three critical points 
-{2N - I)c2fi2/N, 0, and (27V - I)c2fi2/{N + 1) (corre- 
sponding to —4, 0, and 4 in Fig. [1]). The extreme states 
\Fi,F2,F,F) are classified into FF, MM_, MM+, and 
AA phases as in Table 1. Other degenerate states are 
found by repeated applications of Fi^ + -F2-, 

\Fi,F2,F,m) = {Fi^ + F2-f-"'\Fi,F2,F,F) , (7) 

with m = 0,±1,---±F. 

The ground state for a spin-1 polar condensate is frag- 
mented ^fl] , described by a spin singlet state of the form 
(^t)JV/2 |Q^_ ^t (^) 13 invariant under rotations and 
commutes with F2 and F2z- For ferromagnetic conden- 
sate, the condensate ground state favors all atoms aligned 
along the same direction, i.e., takes the form {a[)^ |0), 
and is more stable. When mixing the two together, we 
expect the polar atoms are strongly affected, while the 
back action on to the more stable ferromagnetic atoms is 
negligible. 

The FF phase may be simply described as 
Z^''^{a\)^{h\)'^\Q), with ah polar atoms slaved 
into the same direction as the ferromagnetic ones. 
For -(2A^ - \)c2l32/N < ci2/3 < 0, the MM_ phase 



arises when polar atoms are partly polarized in the 
same direction with the ferromagnetic ones, as if 
there were a finite B-field. Increase of the cou- 
pling interaction (|ci2/3|) breaks singlet pairs in polar 
atoms one by one with AF2 — 2, and gives rise 
to stepwise increases of the total spin. Similar to 
the FF phase, the MM_ phase has a simpler form, 
Zi/2(at)W(5t)i=^2(st)(A^-f2)/2|o)^ with its energy 
E = ciPiNiN + l)/2 + C2(32F2iF2 + l)/2 + C12PNF2/2. 
This phase corresponds to the exact grou nd state for 
F2 - 1< -ci2/3iV/2c2/32 - 1/2 < F2 + 1 III, and we find 
the order parameters 



{Fl) - 



{Fi) 



N{N + l), 
N ^' 



(Fl • F2 



2C2/32 

N^ 

2^2 



(C12/3) 

C12/3 



1 

4' 



N 



(8) 



agree with the mean-field results [Ijl to terms ~ 1/A^. 

For < C12/3 < (27V- l)c2/32/(7V + 1), however, the 
MM^ phase favors polar atoms polarized opposite to the 
ferromagnetic atoms resulting in a decreased total spin. 
This situation is more complicated because all states sat- 



isfying mi +m2 = N — F2 are involved. The ground state 
energy is E = ci(3iN{N + l)/2 + 02/33^^2(^^2 + l)/2 - 
ci2/3(iV + l)F2/2, and we find tliat 



We find interestingly that the total spin vanishes, while 
the species-spins satisfy {F-}) = {F2) ~ N{N + 1) (see 
Fig. [T]). 



(Ff) ^ N{N + 1), 



{Fi) 



{Fi-F2) = 



f N + 1 

V 2C2/32 



(C12/3)' 



{N + lf ^ N+1 
^^^^ci2/3-^-. 



(9) 



The stepwise fine structure of the order parameters is 
not included in Eqs. ^ and ^. When Ci2f3 > 
{2N — 1)02^2/ [N + 1), again all states satisfying mi + 
?7i2 = are included and the AA phase is described by a 
singlet state |18j . 



|iV,iV,0,0) 



N 

E 



IV. ATOMIC NUMBER FLUCTUATIONS 

As a special case, the AA phase is a singlet state, 
which enables us to derive analytically the average par- 
ticle numbers and their associated fluctuations. With 
the help of the methods developed in a related cav- 
ity QED problem 24], we expand the eigenvectors 
I A^, mi) in terms of the Fock states which are defined as 



n\i^ 



n ,4'\"^]) 



,0) 



,0) 



„(i) „(i) „(i)\ 

"■1 : "-0 ^ "--1 / 



a 



0,±1. 



^N m -N -m 1^' ''^i) 1"^' ~''^i) ■ (10) Explicitly, for even and odd mi, we find respectively 





(W-|™i|)/2 


|A^,mi) = 


2^ ^rmk 




fc=0 




(A'-|mi|-l)/2 


|A^,mi) = 


/ J n 



^±^.k,2k,^^^^~k 



m\k 



fc^O 



^ + "-^"^-fc,2fc-H,^-"-^-^-fc 



(11) 



r 



where. 



B 



mik 



[(2fc)!(i^^ - fc)!(^^JjiiJ- - fc)!]i/2' 



(,/7i\2k+l( (N+mi)\(N-mi)l. ^l/2 n\n 
_gV"; ^ v^^^ V 27V! > ^^^ ■ (^2) 



We calculate the particle numbers and number fiuctua- 
tions in the AA phase and find that the average numbers 
of atoms in the six components are exactly all equal, 



4^-)) 



M) 



.^]) 



A^/3, 



(13) 



i.e., the condensate is fragmented [iO|. The fiuctuations 
are given explicitly 



An, 



U) 



VWT9N 
3V5 ' 



U)\ 2^N^ + 3N/2 



(14) 



which approximatively satisfy (An^j = 2/AnQ \ = 

An_{) for large N, as opposed to 2 (Ani) = (Ano) = 
2 (An_i) for the single species singlet state (c2 > 0) [W^ . 



For a comprehensive understanding of the fluctua- 
tions in the entire parameter region, we consider the 
direct product of the Fock states of the two species 

, which may be equiv- 



(1) (1) (1) 



alently defined as 



(2) (2) (2) 

i\ ',n^ ',n^i 



l(l,2) 



(1) (2) 

n^ ',mi,n}^ ' ,m2]m 



(1,2) (1) (2) 

n^' '' n^j %mi,nj, ',m2;m 



(15) 



Here m,i and m2 are the corresponding magnetization 
specified as rrij 



Aj) 



Aj) 



and m 



mi 



7712 IS 



the total magnetization. For simplification, we restrict 
ourselves into the subspace that the total magnetization 
is conserved 771 — 0, in which case all states are non- 
degenerate. Using the full quantum approach of exact 
diagonalization, we simulate the distribution in this sub- 
space with iVi = iV2 — 100, and illustrate the dependence 
of the particle numbers and fluctuations for the six com- 
ponents on C12/3 in Fig. [2] 

For very large values of |ci2/3|, the distribution of 
atoms and fluctuations behave uniformly, indicating 
typical ferromagnetic(FF)/anti-ferromagnetic(AA) spin- 
exchange features. The number distributions for all com- 
ponents are essentially the same over the entire region, 
except for the case of Ci2fi = 0. A tiny C12P ^ 
brings the system into the ferromagnetic-like distribu- 




FIG. 2: (Color online) The dependence of atom numbers and 
fluctuations on ci2/3 at fixed values of ci/3i = — f and 02/^2 = 

2. Blue solid (red dashed) lines denote (na) {[Aria ))■ 

Note /n^^A = /nL'l) and (An'fl\ = /AnL'A. Left (right) 
column denotes the ferromagnetic (polar) condensate. 



,W)\ - /JJ)\ 



,W\ - 



con- 



tion, i.e. 2(nY') = (n'^") = 2(nL'() = N/2, 

sistent with the prediction that the ferromagnetic con- 
densate is more stable. The anti-ferromagnetic distri- 



= (n']l{), in agreement with the analytical re- 



butions for the six components are the same 

^U)\ _ /Jj)\ 

suits (HH). The fluctuations are, on the other hand, 
quite different for positive or negative ci2/3. In the FF 
phase (ci2/3 < —4) the fluctuations for both species 
are small (~ VN). In the AA phase (ci2/3 > 4) 
the fluctuations are large and approximatively satisfy 



An{'^) = 2 / Ani''>) = ( An 



'1 / 



and ( An, 



Af)\ 



{with (An^J.]) « 30.03 



, ~ 15.56). In the regions of MM_ and MM+ 

(—4 < C12/3 < 4), polar atoms exhibit larger fluctuations 
than the ferromagnetic ones, which change quadratically 
with C12/3. When ci2/3 > 0, as the inter-species cou- 
pling increases, the fluctuation of polar atoms dramati- 



cally changes from 2 ( An\ ' ) = (A' 



,(2) 



.(2) 



2(An^^l)to 



An 



(2) 



2 (An, 



.(2) 



An_[) (approximatively). 



Before concluding, we show the difference of the 
singlet state |A^, TV, 0, 0) and the fully-paired state 
Z^'^{&\r,) ^ 1 0) from the viewpoint of quantum fluctu- 
ation [13, [2QI . Using the method of generating function 
[18| . the difference between the two states can be shown 
by taking A^i = A^2 = 2 as an example. The result re- 



mains the same from the angular momentum theory, i.e. 
12,2,0,0)= E <1^'i'?F.%|2,™i)|2,m2) 



^^iell-lA^Brm 



(16) 



More generally, according to the multinomial theorem 

n k 

{xi +X2 + X3)" = Xm c„/fea;5'"'=a;^~'4' (17) 
k=a 1=0 

with Cnik = nl/ {ll{k — l)\{n — k)\), we find that the 
state (012)^ |0) can be described by the Fock state 



(1) (1) (1) 
n\ ,n}) ,n'_{ 



(2) (2) (2)\ 

ni ,ng ,n_i ) as 



{QLr io) 



N k 



^tf;t^w-fe/ ;;tj;t \k-if_;,f luii 



k=Q 1=0 

N k 

= E E^^^)*"^' \k-l,N-k,l)(»\l,N-k,k- I) ,(18) 



k=o 1=0 



, w 



where we have used the property (d^) |0) — ^/n^.\N). 
We calculate the atom numbers and fluctuations for the 
two states |iV, A^,0,0) and Z^/'^{e\2)'^ |0), and find that 
one cannot distinguish them from the atom number dis- 
tributions 

(n(^")) = (n(^")) = (n(^])=iV/3, (19) 

which are the same. The fiuctuations, however, reveal the 
secret. For state Z^/^(6{2)^ |0)i the number fluctuations 
are equally distributed, i.e. 



An 



0)\ 



An 



(j) 



An 



0)^ 



y/N{N + l)/6-N^/9, (20) 



different from the results for the state | A^, A^, 0, 0) we ob- 
tained in Eq. P^ . 



V. CONCLUSION 

To conclude, we study atom number distributions and 
fluctuations for the ground state of a mixture of two spin- 
1 condensates, one being ferromagnetic and the other 
being polar, in the absence of a B-field. For all possible 
inter-species coupling parameter Ci2/9, the exact quantum 
states are constructed from angular momentum algebra 
which are further expanded into the six-component Fock 
states of the mixture. The ground state quantum phases 
are classifled into four types according to ci2/3. The most 



interesting AA phase of a singlet for the total spin, where 
spins of each species are polarized in opposite directions, 
is fragmented with six components equally distributed 
(n-f) — (uq ) = \n_i) = N/3 and exhibits anoma- 
lous number fluctuations. We find that the fluctuations 



satisfy (An^^'^) = 2 ( Anf^'^ 



An^J\ for large N, dif- 
ferent from the result of single species polar state (c2 > 0) 
2 (Ani) = (Ano) = 2 (An_i). Our results highlight the 
significant promises for experimental work on Na and Rb 
atomic condensate mixtures since optical Feshbach reso- 



nances make it possible to tune the spin exchange inter- 
action. 
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